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Abstract. I consider novel ways by which Gamma- Ray Bursts (GRBs) and quasars 
can be used to probe the universe. Clues about how and when was the intergalactic 
medium ionized can be read off the UV emission spectrum ol GRB explosions Irom the 
first generation ol stars. The existence ol intergalactic and galactic stars can be interred 
from their gravitational microlensing effect on GRB afterglows. Prior to reionization, 
quasars should be surrounded by a halo ol scattered Lya radiation which probes the 
neutral intergalactic medium (IGM) around them. The situation is analogous to the 
appearance ol a halo ol scattered light around a street lamp which is embedded in a 
dense log. Outflows from quasars magnetize the IGM at all redshitts. As a result, the 
shocks produced by converging flows during the formation ol large scale structure in the 
IGM, accelerate electrons to relativistic energies and become visible in the radio regime 
through their synchrotron emission and in the 7-ray regime through their inverse- 
Compton scattering ol the microwave background photons. During transient episodes 
ol strong mergers, X-ray clusters should therefore appear as extended radio or 7-ray 
sources on the sky. 



1 Novel Cosmological Studies with Gamma-Ray Bursts 
1.1 Preface 

Since their discovery four decades ago, quasars have been used as powerful light- 
houses which probe the intervening universe out to high redshifts, z ~ 6 ]4^pt| . 
The spectra of almost all quasars show strong emission lines of metals, indicating 
super-solar enrichment of the emitting gas [M . This implies that at least in the 
cores of galaxies, formation of massive stars and their evolution to supernovae 
preceded the observed quasar activity. If Gamma-Ray Bursts (GRBs) originate 
from the remnants of massive stars (such as neutron stars or black holes), as 
seems likely based on recent estimates of their energy output |89| , p4 123 1 , then 



they should exist at least out to the same redshift as quasars. Although GRBs are 
transient events, their peak optical-UV flux can be as bright as that of quasars. 
Hence, GRBs promise to be as useful as quasars in probing the high-redshift 
universe. 

Not much is known observationally about the universe in the redshift interval 
z = 6-30, when the first generation of galaxies condensed out of the primordial 
gas left over from the Big Bang (see reviews [|| and |36|). Observations of the 
microwave background anisotropics indicate that the cosmic gas became neutral 
at z ~ 1000 and remained so at least down to z ~ 30 (see, e.g. j87[). On the 
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other hand, the existence of transmitted flux shortward of the Lya resonance 
in the spectrum of the highest-redshift quasars and galaxies (see, for example, 
Figure , indicates that the intergalactic medium was reionized to a level better 
than 99.9999% by a redshift z ~ 6. This follows from the fact that the Lya 
optical depth of the intergalactic medium at high-redshifts (z 3> 1) is pSj ], 

Tre 2 f a X a n HI (z) 5 / Q b h\ ( D m \ ~ 1/2 ( 1 + z x 3 ■' 
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where i? « lOO/i km s" 1 Mpc" 1 J?™ /2 (1 + z) 3 / 2 is the Hubble parameter at the 
source redshift z, f a = 0.4162 and A Q = 1216A are the oscillator strength and 
the wavelength of the Lya transition; tihi(z) is the average intergalactic density 
of neutral hydrogen at the source redshift (assuming primordial abundances); 
fi m and fib are the present-day density parameters of all matter and of baryons, 
respectively; and xhi is the average fraction of neutral hydrogen. Modeling |^o| 
of the transmitted flux in Figure ^ implies r Q < 0.5 or xhi < 10~ 6 , i.e., most 
of the low-density gas throughout the universe is ionized to a level of 99.999% 
at z < 6. However, there are some dark intervals in the spectrum which could 
be indicative of regions with a higher neutral fraction In fact, Becker et 
al. Q reported the detection of an extended (> 300A long) dark interval just 
shortward of the Lya emission line in the spectrum of the newly discovered 
quasar SDSS 1030+0524 at z = 6.28. The suppression of the flux by a factor 
> 150 may indicate the first detection of the Gunn-Peterson trough, although 
caution is warranted since the inferred optical depth r Q > 5 can be produced by 
a neutral fraction as small as Xhi ~ 2 x 10~ 4 according to equation ([[]). 

The question: how and when was the universe reionized? defines a new fron- 
tier in observational cosmology J56) . The UV spectrum of GRB afterglows can 
be used to probe the ionization and thermal state of the intergalactic gas during 
the epoch of reionization, at redshifts z ~ 7-10 Q]. The stretching of the tem- 
poral evolution of GRB lightcurves by the cosmological redshift factor (1 + z), 
makes it easier for an observer to react in time and measure a spectrum of their 
optical-UV emission at its peak. 

Energy arguments suggest that reionization resulted from photoionization 
and not from collisional ionization p| , p8| . The corresponding sources of UV 
photons were either stars or quasars. Recent simulations of the first generation 
of stars that formed out of the primordial metal-free gas indicate that these 
stars were likely massive If GRBs result from compact stellar remnants, 

such as black holes or neutron stars, then the fraction of all stars that lead 
to GRBs may have been higher at early cosmic times. This, however, is true 
only if the GRB phenomena is triggered on a time scale much shorter than 
the age of the universe at the corresponding redshift, which for z ^> 1 is ~ 
5.4 x 10 8 yr (V0.7)- 1 (^m/0-3)- 1/2 [(l + z)/10]- 3 / 2 . 

1.2 Properties of High-Redshift GRB Afterglows 

Young (days to weeks old) GRBs outshine their host galaxies in the optical 
regime. In the standard hierarchical model of galaxy formation, the character- 
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Fig. 1. Spectrum of the quasar SDSS1044-0125 at z = 5.73 (Djorgovski et al. 2001), 
originally discovered by the Sloan Digital Sky Survey (Fan et al. 2000). 



istic mass and hence optical luminosity of galaxies and quasars declines with 
increasing redshift |39 79,3|. Hence, GRBs should become easier to observe than 
galaxies or quasars at increasing redshift. Similarly to quasars, GRB afterglows 
possess broad-band spectra which extend into the rest-frame UV and can probe 
the ionization state and metallicity of the IGM out to the epoch when it was 
reionized at redshifts z ~ 7-10 ]5q| . Lamb & Reichart |]54j have extrapolated 
the observed 7-ray and afterglow spectra of known GRBs to high redshifts and 
emphasized the important role that their detection might play in probing the 
IGM. Simple scaling of the long- wavelength spectra and temporal evolution of 
afterglows with redshift implies that at a fixed time lag after the GRB trigger 
in the observer's frame, there is only a mild change in the observed flux at in- 
frared or radio wavelengths as the GRB redshift increases. Ciardi & Loeb [Q 
demonstrated this behavior using a detailed extrapolation of the GRB fireball 
solution into the non-relativistic regime (see the 2/xm curves in Figure |J). De- 
spite the strong increase of the luminosity distance with redshift, the observed 
flux for a given observed age is almost independent of redshift in part because 
of the special spectrum of GRB afterglows (see Figure ^) , but mainly because 
afterglows are brighter at earlier times and a given observed time refers to an 
earlier intrinsic time in the source frame as the source redshift increases. The 
mild dependence of the long- wavelength (A b s > 1/xm) flux on redshift stands 
in contrast to other high-redshift sources such as galaxies or quasars, which 
fade rapidly with increasing redshift j39||7l| |3|| ■ Hence, GRBs provide exceptional 
lighthouses for probing the universe at z = 6-30, at the epoch when the first 
stars had formed. 
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Fig. 2. Theoretical expectation for the observed afterglow flux of a GRB as a function 
of its redshift (from Ciardi & Loeb 2000) . The curves refer to an observed wavelength of 
5000 A(thin lines) and 2fim (thick lines). Different line types refer to different observed 
times after the GRB trigger, namely 1 hour (solid line), 1 day (dotted) and 10 days 
(dashed). The 5000A flux is strongly absorbed at z > 4.5 by intergalactic hydrogen. 
However, at infrared and radio wavelengths the observed afterglow flux shows only a 
mild dependence on the source redshift. 



Assuming that the GRB rate is proportional to the star formation rate and 
that the characteristic energy output of GRBs is ~ 10 52 ergs, Ciardi & Loeb 
| fl2| predicted that there are always ~ 15 GRBs from redshifts z > 5 across the 
sky which are brighter than ~ 100 nJy at an observed wavelength of ~ 2 /mi. 
The infrared spectrum of these sources could be taken with future telescopes 
such as the Next Generation Space Telescope (planned for launch in 2009; see 
http://ngst.gsfc.nasa.gov/), as a follow-up on their early X-ray localization with 
the Swift satellite (planned for launch in 2003; see http://swift.sonoma.edu/). 

The redshifts of GRB afterglows can be estimated photometrically from ei- 
ther the Lyman limit or Lya troughs in their spectra. At low redshifts, the 
question of whether the Lyman limit or Lya trough interpretation applies de- 
pends on the absorption properties of the host galaxy. If the GRB originates 
from within the disk of a star-forming galaxy, then the afterglow spectrum will 
likely show a damped Lya trough. At z > 6 the Lya trough would inevitably 
exist since the intergalactic Lya opacity is > 90% (see Figure 13 in j79j). In- 
terestingly, an absorption feature in the afterglow spectrum which is due to the 
neutral hydrogen within a molecular cloud or the disk of the host galaxy, is likely 
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to be time-dependent due to the ionization caused by the UV illumination of 
the afterglow itself along the line-of-sight jti^j . 

So far, there have been two claims for high-redshift GRBs. Fruchter p5[ 
argued that the photometry of GRB 980329 is consistent with a Lya trough 
due to IGM absorption at z ~ 5. Anderson et al. || inferred a redshift of 
z = 4.5 for GRB 000131 based on a crude optical spectrum that was taken 
by the VLT a few days after the GRB trigger. These cases emphasize the need 
for a coordinated observing program that will alert 10-meter class telescopes to 
take a spectrum of an afterglow about a day after the GRB trigger, based on a 
photometric assessment (obtained with a smaller telecope using the Lyman limit 
or Lya troughs ) that the GRB may have originated at a high redshift. 

In the following two subsections, I illustrate the usefulness of GRB afterglows 
for cosmological studies through two examples. 
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Fig. 3. Sketch of the expected spectrum of a source at a redshift z 3 slightly above 
the reionization redshift z Ic i on . The transmitted flux due to HII bubbles in the pre- 
reionization era, and the Lya forest in the post-reionization era, are exaggerated for 
illustration. 
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1.3 Probing the Reionization Epoch and Metallicity History of the 
IGM 

The UV spectra of GRB afterglows can be used to measure the evolution of 
the neutral intergalactic hydrogen with redshift. Figure || illustrates schemat- 
ically the expected absorption just beyond the reionization redshift. Resonant 
scattering suppresses the spectrum at all wavelengths corresponding to the Lya 
resonance prior to reionization. Since the Lya cross-section is very large, any 
transmitted flux prior to reionization reflects a large volume of ionized hydrogen 
along the line-of-sight. If the GRB is located at a redshift larger by > 18% than 
the reionization redshift, then the Lya and the Ly/3 troughs will overlap. Unlike 
quasars, GRBs do not ionize the IGM around them; their limited energy supply 
~ 10 52 ergs pj|p4|]23| ] can ionize only ~4x 10 5 -M Q of neutral hydrogen within 
their host galaxy. 

Quasar spectra indicate the existence of an IGM metallicity which is a frac- 
tion of a percent of the solar value [JlTj . The metals were likely dispersed into the 
IGM through outflows from galaxies, driven by either supernova or quasar winds 
[||J2£|. Detection of metal absorption lines in the spectrum of GRB afterglows, 
produced either in the IGM or the host galaxy of the GRB, can help unravel 
the evolution of the IGM metallicity with redshift and its link to the evolution 
of galaxies. Detection of X-ray absorption by intergalactic metals can be used to 
establish the existence of the warm component of the IGM which has not been 



observed so far E3,W3 69 



1.4 Cosmological Microlensing of Gamma- Ray Bursts 

Loeb & Perna [|5^] noted the coincidence between the angular size of a solar-mass 
lens at a cosmological distance and the micro-arcsecond size of the image of a 
GRB afterglow. They therefore suggested that microlensing by stars can be used 
to resolve the photospheres of GRB fireballs at cosmological distances. (Alterna- 
tive methods, such as radio scintillations, only provide a constraint on the radio 
afterglow image size [ p4[j8^ | but do not reveal its detailed surface brightness dis- 
tribution, because of uncertainties in the scattering properties of the Galactic 
interstellar medium.) 

The fireball of a GRB afterglow is predicted to appear on the sky as a ring (in 
the optical band) or a disk (at low radio frequencies) which expands laterally at a 
superluminal speed, ~ _Tc, where r ^> 1 is the Lorentz factor of the relativistic 



blast wave which emits the afterglow radiation [p8|j76y68| , |36| . For a spherical 
explosion into a constant density medium (such as the interstellar medium), the 
physical radius of the afterglow image is of order the fireball radius over r, or 
more precisely BfJ 



• j? \ L i° / + \ 5 / 8 
I,'.. :!.!) x !()"•( i-^ j ij+l) <•»>■ l- r > 



where E52 is the hydrodynamic energy output of the GRB explosion in units of 

Jays 



10 52 ergs, «i is the ambient medium density in units of 1 cm 3 , and tdays is the 
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Fig. 4. A typical broken power-law spectrum of a GRB afterglow at a redshift z = 1 
(from Granot & Loeb 2001). The observed flux density, F v , as a function of frequency, 
v, is shown by the boldface solid line at an observed time tdays = 1 for an explosion 
with a total energy output of 10 52 ergs in a uniform interstellar medium (k = 0) with 
a hydrogen density of 1 cm -3 , and post-shock energy fractions in accelerated electrons 
and magnetic field of e c = 0.1 and es = 0.03, respectively. The thin solid line shows the 
same spectrum when it is microlensed by an intervening star with an impact parameter 
equal to the Einstein angle and Ro = [0 S (1 day)/#E] = 1. The insert shows the excess 
evolution of the break frequencies v\> = v a , v m and v c (normalized by their unlensed 
values) due to microlensing. 

observed time in days. At a cosmological redshift z, this radius of the GRB image 
occupies an angle 9 S = R s /Da, where Da(z) is the angular diameter distance 
at the GRB redshift, z. For the typical cosmological distance, Da ~ 10 28 cm, 
the angular size is of order a micro-arcsecond (/ias). Coincidentally, this image 
size is comparable to the Einstein angle of a solar mass lens at a cosmological 
distance, 

where Mi ons is the lens mass, and D = D 0S D \/ D\ s is the ratio of the angular- 
diameter distances between the observer and the source, the observer and the 
lens, and the lens and the source (77). Loeb & Perna J5?J showed that because 
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the ring expands laterally faster than the speed of light, the duration of the 
microlensing event is only a few days rather than tens of years, as is the case 
for more typical astrophysical sources which move at a few hundred km s" 1 or 
- lQ- 3 c. 

The microlensing lightcurve goes through three phases: (i) constant magni- 
fication at early times, when the source is much smaller than the source-lens 
angular separation; (ii) peak magnification when the ring-like image of the GRB 
first intersects the lens center on the sky; and (iii) fading magnification as the 
source expands to larger radii. 




Fig. 5. The surface brightness, normalized by its average value, as a function of the 
normalized radius, r, from the center of a GRB afterglow image (where r = at the 
center and r — 1 at the outer edge). The image profile changes considerably between 
different power-law segments of the afterglow spectrum, F v oc v 13 (see Figure [ij). There 
is also a strong dependence on the power-law index of the radial density profile of the 
external medium around the source, p oc R' k (taken from Granot & Loeb 2001). 
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Granot & Loeb |35| calculated the radial surface brightness profile (SBP) of 
the image of a Gamma-Ray-Burst (GRB) afterglow as a function of frequency 
and ambient medium properties, and inferred the corresponding magnification 
lightcurves due to microlensing by an intervening star. The afterglow spectrum 
consists of several power-law segments separated by breaks, as illustrated by Fig- 
ure § The image profile changes considerably across each of the spectral breaks, 
as shown in Figure 0. It also depends on the power-law index, k, of the radial 
density profile of the ambient medium into which the GRB fireball propagates. 
Gaudi & Loeb |nj have shown that intensive monitoring of a microlensed af- 
terglow lightcurve can be used to reconstruct the parameters of the fireball and 
its environment. The dependence of the afterglow image on frequency offers a 
fingerprint that can be used to identify a microlensing event and distinguish it 
from alternative interpretations. It can also be used to constrain the relativistic 
dynamics of the fireball and the properties of its gaseous environment. At the 
highest frequencies, the divergence of the surface brightness near the edge of the 
afterglow image (r = 1 in Figure ||) depends on the thickness of the emitting 
layer behind the relativistic shock front, which is affected by the length scale 
required for particle acceleration and magnetic field amplification behind the 
shock[|60|]37il. 



Ioka & Nakamura 44 considered the more complicated case where the ex- 
plosion is collimated and centered around the viewing axis. More general orien- 
tations that violate circular symmetry need to be considered in the future. 



GRB 000301C Garnavich, Loeb, & Stanek ]2^] have reported the possible 
detection of a microlensing magnification feature in the optical-infrared light 
curve of GRB 000301C (see Figure ||). The achromatic transient feature is well 
fitted by a microlensing event of a 0.5M@ lens separated by an Einstein angle 
from the source center, and resembles the prediction of Loeb & Perna 57 for a 
ring- like source image with a narrow fractional width (~ 10%). Alternative inter- 
pretations relate the transient achromatic brightening to a higher density clump 
into which the fireball propagates ||, or to a refreshment of the decelerating 
shock either by a shell which catches up with it from behind or by continu- 
ous energy injection from the source [ jjjl) . However, the microlensing model has 
a smaller number of free parameters. If with better data, a future event will 
show the generic temporal and spectral characteristics of a microlensing event, 
then these alternative interpretations will be much less viable. A galaxy 2" from 
GRB 000301C might be the host of the stellar lens, but current data provides 
only an upper-limit on its surface brightness at the GRB position. The existence 
of an intervening galaxy increases the probability for microlensing over that of 
a random line-of-sight. 

Gaudi, Granot, & Loeb |52| have shown that direct inversion of the observed 
light curve for GRB 000301C yields a surface brightness profile (SBP) of the 
afterglow image which is strongly limb-brightened, as expected theoretically (see 
Figure 0) . 
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Fig. 6. UBVRIJK photometry of GRB 000301C as a function of time in days from 
the GRB trigger (from Garnavich et al. 2000; Gaudi et al. 2001). Data points have 
been offset by the indicated amount for clarity. The dashed line is the best-fit smooth, 
double power- law lightcurve (with no lensing) , while the solid line is the overall best-fit 
microlensing model, where the SBP has been determined from direct inversion. The 
inset shows the i?-band data only. The dotted line is the best-fit microlensing model 
with theoretically calculated SBP, for k — and v > v c - 

Obviously, realistic lens systems could be more complicated due to the exter- 
nal shear of the host galaxy or a binary companion. Mao & Loeb |53| calculated 
the magnification light curves in these cases, and found that binary lenses may 
produce multiple peaks of magnification. They also demonstrated that all after- 
glows are likely to show variability at the level of a few percent about a year 
following the explosion, due to stars which are separated by tens of Einstein 
angles from the line-of-sight. 

What is the probability for microlensing? If the lenses are not strongly clus- 
tered so that their cross-sections overlap on the sky, then the probability for 
having an intervening lens star at a projected angular separation 9 from a source 
at a redshift z ~ 2 is - O.3J2*(0/0 B ) 2 0,§|56],§7), where f2* is the cosmological 
density parameter of stars. The value of .!?* is bounded between the density of 
the luminous stars in galaxies and the total baryonic density as inferred from 
Big Bang nucleosynthesis, 7 x 10~ 3 < < 5 x 1CT 2 j27|. Hence, all GRB af- 



Novel Ways to Probe the Universe with Gamma-Ray Bursts and Quasars 



11 



w 

CD 

si: 
I 

CD 

o 

< 

: : 

m 

c; 
c:.; 
x 

a 

o 



5 
4 
3 
2 
1 



5 

4 
3 
2 

: 

o 



-(a) 



Mpn Direct Inversion, k = 0, ^ a =149 (1(7 errors) j- 
{3o" errors) 

k = 0, V.<V<V m , x 2 =190 

k = 0, v m <v<v v x 2 =166 

k = 0, v>v„ x 2 =155 



/ 





/ 



■ 



-^r — I 



L (b) 



h-^H Direct Inversion, k = S, ^ s =147 (1(7 errors) 

(3(7 errors) -i — f- 

k = 2, i/ i <iy<y ni , ^=217 | -i- 

k = 2, i/ m <i/<i/., x 2=19a 

k = 2, 1/>K , x z =156 



7 
./ 



/ 



"■■■\ 



0.5 

Normalized Radius 



Fig. 7. Fitting GRB 000301C with different SBPs as a function of normalized radius 
(taken from Gaudi et al. 2001). The points are the SBPs determined from direct in- 
version, with la errors (solid) and 3cr errors (dotted). The curves are theoretically 
calculated SBPs for various frequency regimes (see Figure ||). (a) Uniform external 
medium, k = 0. (b) Stellar wind environment, k — 2. The number of degrees of free- 
dom is 92 for the direct inversion points and 89 for the curves. 



terglows should show evidence for events with 9 ~ 300 e, for which microlensing 
provides a small perturbation to the light curve |33|. (This crude estimate ig- 
nores the need to subtract those stars which are embedded in the dense central 
regions of galaxies, where macrolensing dominates and the microlensing optical 
depth is of order unity.) However, only one out of roughly a hundred afterglows 
is expected to be strongly microlensed with an impact parameter smaller than 
the Einstein angle. Indeed, Koopmans & Wambsganss Esl have found that the 
'a posteriori' probability for the observed microlensing event of GRB 000301C 
along a random line-of-sight is between 0.7-2.7% if 20-100% of the dark matter 
is in compact objects. 

Microlensing events are rare but precious. Detailed monitoring of a few strong 
microlensing events among the hundreds of afterglows detected per year by the 
forthcoming Swift satellite, could be used to constrain the environment and the 
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dynamics of relativistic GRB fireballs, as well as their magnetic structure and 
particle acceleration process. 

2 Illumination of the Intergalactic Medium by Quasars 

2.1 Lya Halos 

The absorption trough created in the spectrum of a distant sources due to neu- 
tral hydrogen diminishes down to undetectable flux levels as soon as the mean 
neutral fraction of hydrogen is larger than ~ 10~ 4 (see equation |l|). Once the 
transmitted flux reaches very low levels, it is difficult to infer whether the neutral 
fraction of hydrogen is as low as 0.01% or as high as 100%. A novel way to study 
the reionization phase transition relies on a direct detection of intergalactic neu- 
tral hydrogen. Loeb & Rybicki |38|,[75) have shown that the neutral gas prior 
to reionization can be probed through narrow-band imaging of embedded Lya 
sources. The physical situation is analogous to the appearance of a halo of scat- 
tered light around a street lamp which is embedded in a dense fog. The "street 
lamp" in this metaphor is a high-redshift quasar which emits Lya photons into 
the surrounding neutral hydrogen gas. The IGM scatters these photons and acts 
as a fog. 

The radiation of the first galaxies is strongly absorbed shortward of their rest- 
frame Lya wavelength by neutral hydrogen in the intervening IGM. However, 
the Lya photons emitted by these sources are not eliminated but rather scatter 
until they redshift out of resonance and escape due to the Hubble expansion 
of the surrounding intergalactic hydrogen (see Fig. 8). Detection of the diffuse 
Lya halos around high redshift sources would provide a unique tool for prob- 
ing the neutral IGM before the epoch of reionization. Loeb & Rybicki ]5§| , |75| 
explored the above effect for a uniform, fully-neutral IGM in a pure Hubble 
flow. It is important to extend this analysis to more realistic cases of sources 
embedded in an inhomogeneous IGM, which is partially ionized by these sources 
[Eo|. It would be interesting to extract particular realizations of the perturbed 
IGM around massive galaxies from numerical simulations, and to apply a suit- 
able radiative transfer code in propagating the Lya photons from the embedded 
galaxies. Mapping of the properties of the associated Lya halos will allow to 
assess their detectability with future observations. 

2.2 Magnetization of the Intergalactic Medium by Quasar Outflows 

Magnetic fields and cosmic rays play an important dynamical role in the in- 
terstellar medium of galaxies. However, they are often ignored in discussions of 
the intergalactic medium (IGM). How significant is the pressure from these non- 
thermal components in the IGM? Has the magnetic field observed in collapsed 
objects, such as galaxies or clusters of galaxies, originated from the IGM? The 
first question is particularly relevant for the reconstruction of the mass distribu- 
tion of galaxy clusters from X-ray data. 
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Fig. 8. Loeb-Rybicki halos: Scattering of Lya line photons from a galaxy embedded in 
the neutral intergalactic medium prior to reionization. The line photons diffuse in fre- 
quency due to the Hubble expansion of the surrounding medium and eventually redshift 
out of resonance and escape to infinity. A distant observer sees a Lya halo surrounding 
the source, along with a characteristically asymmetric line profile. The observed line 
should be broadened and redshifted by about one thousand km s _1 relative to other 
lines (such as H a ) emitted by the galaxy. 



Outflows from quasars inevitably pollute the IGM with magnetic fields. The 
short-lived activity of a quasar leaves behind an expanding magnetized bubble in 
the IGM. Furlanetto & Loeb |^8| modeled the expansion of the remnant quasar 
bubbles and calculated their distribution as a function of size and magnetic field 
strength at different redshifts. They have found that generically by a redshift 
z ~ 3, about 5-20% of the IGM volume is filled by magnetic fields with an 
energy density > 10% of the mean thermal energy density of a photo-ionized 
IGM at ~ 10 4 K (see Figure ||). As massive galaxies and X-ray clusters condense 
out of the magnetized IGM, the adiabatic compression of the magnetic field 
could result in the field strength observed in these systems without a need for 
further dynamo amplification. The intergalactic magnetic field could also provide 
a non-thermal contribution to the pressure of the photo-ionized gas that may 
account for the claimed discrepancy between the simulated and observed Doppler 
width distributions of the Lya forest at z > 2 |l^,Q. However, the supplied 
magnetic energy is unlikely to be dynamically important today since the present- 
day IGM was heated by gravitationally-induced shocks to an average mass- 



weighted temperature of ~ 3 x 10 K 13 , larger by two orders of magnitude 



than the magnetic energy input from quasars. 
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Fig. 9. Upper Panel: Volume filling fraction of magnetized bubbles F(z), as a function 
of redshift (taken from Furlanetto & Loeb 2001). Lower Panel: Ratio of magnetic 
energy density, us/e-i, to a fiducial thermal energy density Ufid = 3n(z)kTiGM of a 
photoionized IGM, where Tigm = 10 4 K, as a function of redshift. The solid curves 
assume a different minimum mass for a galaxy out of which a quasar outflow may 
originate, Mh,min- Examples include a case where this minimum mass is determined 
by atomic cooling before reionization and by infall suppression afterward (top curve), 
and fixed-value cases with Mh,min = 1O 9 M (middle curve), and Mh,min = 1O 1O M0 
(bottom curve). The vertical dotted line indicates the assumed redshift of reionization, 

Z r = 7. 

2.3 Particle Acceleration in Magnetized Intergalactic Shocks 

Even though the present-day magnetic pressure may not be dynamically im- 
portant, the magnetization of the IGM by quasar outflows has important conse- 
quences. Once magnetized, the shocks produced in the IGM during the formation 
of large-scale structure can accelerate a population of highly relativistic electrons 
by the Fermi mechanism, similarly to the collisionlcss shocks of supernova rem- 
nants. The accelerated electrons emit synchrotron radiation as they gyrate in 
the embedded magnetic fields and produce 7-ray radiation by inverse-Compton 
scattering the microwave background photons. This non-thermal radiation spans 
some ~ 20 orders of magnitude in photon energies from the radio to the TeV 
7-ray regime. The radiation is expected to delineate the strong shocks formed 
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in the cosmic web of large scale sheets and filaments in the IGM. The bright- 
est emission originates from the strongest intergalactic shocks around galaxy 
clusters. 

2.4 Non-Thermal Radio and Gamma-Ray Emission 

More than a third of all X-ray clusters with luminosities > 10 45 erg s _1 pos- 
sess diffuse radio halos Q. Based on energy arguments and circumstantial ev- 
idence, these radio halos are believed to be caused by synchrotron emission 
from shock-accelerated electrons during the merger events of their host clusters 
( |4l] , ^ , pl| ; see [||| for references to alternative, less successful models). These 
highly-relativistic electrons cool primarily through inverse- Compton scattering 
off the microwave background. Since their cooling time is much shorter than the 
dynamical time of their host cluster, the radio emission is expected to last as 
long as the shock persists and continues to accelerate fresh electrons to rela- 
tivistic energies. Intergalactic shocks also occur along the filaments and sheets 
that channel mass into the clusters. These structures, also traced by the galaxy 
distribution Jl6| , are induced by gravity and form due to converging large-scale 
flows in the IGM. 

The intergalactic synchrotron emission contaminates the cosmic microwave 
background (CMB) anisotropies - relic from the epoch of recombination, and as 
such needs to be considered in the design and analysis of anisotropy experiments 
at low frequencies. Previous estimates of synchrotron contamination of the cos- 
mic anisotropies focused on Galactic emission, which occurs primarily at low 
Galactic latitudes and large angular scales Q . Waxman & Loeb Q calculated 
the intergalactic synchrotron contribution to the fluctuations in the radio sky as 
a function of frequency and angular scale (see Figure |l^). They assumed that 
most of the emission originates from the virialization shocks around X-ray clus- 
ters, and used the Press-Schechter Q mass function to describe the abundance 
of such clusters as a function of redshift (see also subsequent work in Q). Al- 
though the synchrotron background amounts to only a small fraction of the CMB 
intensity, Waxman & Loeb [|90f showed that its fluctuations could dominate over 
the primordial CMB fluctuations at low photon frequencies, v < 10 GHz. They 
found that radio emission from cluster shocks contributes a fluctuation ampli- 
tude of ~ 40/iK x(^/10GHz)~ 3 to the CMB on angular scales between 1 and 
0.1°, respectively. Interestingly, current anisotropy experiments are just sensitive 
to this level of fluctuations. Existing detections by CAT (50 ± 15/j.K at 15 GHz 
on 0.2-0.5° scales) and OVRO (56lJ;^K at 20 GHz on 0.1-0.6° scales), as well 
as 95% upper limits (< 40/iK on arcminute scales at 9-15 GHz by the ATCA, 
RYLE and VLA detectors) are consistent with our prediction. 

Loeb & Waxman |39| have shown that a significant fraction of the diffuse 
7-ray background |7§| ] might have been generated by the shocks resulting from 
the formation of large-scale structure in the intergalactic medium. Similarly to 
the collisionless shocks of supernova remnants , these shocks produce a popu- 
lation of highly-relativistic electrons with a maximum Lorentz factor > 10 7 that 
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Fig. 10. Fractional intensity fluctuations, {{I V {Q)I V {6)) — (ij,) 2 ) 1 / 2 /(/„), in synchrotron 
(radio) and inverse-Compton (7-ray) background flux (taken from Waxman & Loeb 
2000). The dimensionless coefficient, / s h, is of order unity. The ratio of synchrotron 
to CMB intensity is (Il yn )/I? MB = 6 x 10" 6 / s h 2 (?s/0.01)(i//10GHz)- 3 , where the 
post-shock magnetic energy fraction £b is assumed to be of order 0.01, corresponding 
to a magnetic field strength of 0.1/xG in the virialization shock of X-ray clusters. The 
post-shock energy fraction carried by relativistic electrons is assumed to be £ e ~ 0.05, 
similar to the value inferred in the shocks of supernova remnants. 

scatter a small fraction of the microwave background photons in the present- 
day universe up to 7-ray energies, thereby providing the 7-ray background. The 
predicted flux agrees with the observed diffuse background over more than four 
decades in photon energy and is not sensitive to the precise magnetic field value, 
provided that the fraction of shock energy carried by relativistic electrons is 
£ e ~ 0.05 (a value consistent with that inferred for supernovae shocks). The 
same electrons that emit 7-rays by inverse-Compton scattering of microwave 
background photons, also produce synchrotron radiation in the radio band due 
to intergalactic magnetic fields. The existence of magnetic fields with an ampli- 
tude > 0.1/iG, is inferred in cluster halos p^ , |30]|74||45| ] , and is also required for 
the Fermi acceleration of these electrons. The appearance of radio halos around 
young X-ray clusters is therefore a natural consequence and an important test 
of this model for the extragalactic 7-ray background. The combination of radio 
and 7-ray data can be used to calibrate £ e and determine the strength of the 
intergalactic magnetic field. 

The assumed acceleration of electrons by collisionless shocks is similar to 
that observed in supernova remnants. Recent X-ray [^9|j50t| and TeV J8(|^35) ob- 
servations of the supernova remnants SN1006 and SNR RX J1713. 7-3946 imply 
that electrons are accelerated in the remnant shocks up to an energy ~ 100 TeV, 
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and are confined to the collisionless fluid by magnetic fields. These shocks have 
a velocity of order 10 3 km s _1 , similar to the velocity of the intergalactic shocks 
we consider here. Although the plasma density is very different in the two cases, 
the density may be scaled out of the problem by measuring time in units of 
ui~^ , where cu pe is the electron plasma frequency (shock characteristics may also 
depend on the pre-shock magnetic field, which introduces a dimensionless pa- 
rameter into the problem, uj ce /u pe , where uj cb is the electron cyclotron frequency. 
However, for both supernova and IGM shocks, uj ce /uj pe <C 1, suggesting a simi- 
lar behavior in both cases). For supernova shocks, the inferred energy density in 
relativistic electrons constitutes 1-10% of the post-shock energy density in these 
remnants fisfl , a fraction consistent with the global ratio between the mean en- 
ergy density of cosmic-ray electrons and the turbulent energy density in the 
interstellar medium of our galaxy. 

The production of anisotropic backgrounds of radio and 7-ray radiation by 
strong intergalactic shocks is a natural consequence of structure formation in 
the Universe. The brightest emission originates from the shocks on Mpc scales 
around newly formed, massive X-ray clusters. The foreground synchrotron fluc- 
tuations might be comparable to the anisotropy signals detected by existing 
low-frequency microwave background experiments, and can be easily isolated 
through multi-frequency observations. Polarization anisotropy experiments could 
then constrain the coherence length of the intergalactic magnetic field. 

Different astrophysical sources may contribute to the strength of the inter- 
galactic magnetic fields; aside from radio sources |l5],^8| they include: the large 
scale shocks themselves J52|, the first generation of stars |73| and supernova- 
driven winds from galaxies [Q. Different contributions may lead to different 
scalings of the magnetic field energy density with IGM gas parameters. The non- 
thermal 7-ray emission is only weakly dependent on the magnetic field strength 
[ |59| . However, the non-thermal radio emission is highly sensitive to this parame- 
ter and different models may lead to a different dependence of the radio emission 
on scale and redshift, which may have observable consequences (e.g. modifying 
the functional dependence of the correlation function shown in Figure |Io|) . 

Recent observations of the giant radio galaxy NGC 315 may constitute the 
first direct detection of collisionless large scale structure shocks, accelerating elec- 
trons to high energies jl9| . Radio observations of this source are most naturally 
explained by the existence of a large scale shock produced by a flow converging 
towards an intersection of galaxy filaments |l9|| . 

Waxman & Loeb |)(J predicted a fluctuation amplitude > 40% in the 7-ray 
background intensity on sub-degree scale, and the existence of extended, > 1°, 7- 
ray halos, associated with newly formed massive clusters. On scales larger than 
a degree the fluctuation amplitude declines and is well below the anisotropy 
limits from EGRET (see Figure 5 in (78)). Detection of the predicted signals 
will provide a calibration of the uncertain model parameter, £ e . The high-energy 
maps required to detect the predicted anisotropy signal will be made between 
20 MeV and 300 GeV by the G LAST instrument (planned for launch in 2005; 
see http://glast.gsfc.nasa.gov/), which is expected to be more sensitive than 
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EGRET by an order-of-magnitude Q . The predicted 7-ray halos may constitute 
a significant fraction of the unidentified extra-Galactic EGRET sources p^j85[ |. 
However, since the angular extension of the brightest halos is large, a more 
careful analysis is required in order to assess the detectability of such halos by 
EGRET. 

A future, dedicated, all-sky anisotropy experiment, operating at several fre- 
quencies below 10 GHz, would be able to map the fluctuations in the inter- 
galactic synchrotron background. The resulting synchrotron map could then be 
cross-correlated with full-sky maps at hard X-ray or 7-ray energies to confirm its 
cosmic origin. Identification of the synchrotron fluctuations together with their 
counterpart inverse-Compton emission of hard X-rays or 7-rays by the same pop- 
ulation of shock-accelerated electrons, can be used to empirically determine the 
strength and spatial distribution of the intergalactic magnetic field. Similarly, 
the correlation between radio and 7-ray halos may be detectable around indi- 
vidual X-ray clusters. Strong radio halos could be the best indicators for bright 
7-ray clusters, which would provide the first obvious targets for GLAST. 

Semi-analytic models ]9C| identify the intergalactic shocks with smooth spher- 
ical accretion of gas onto clusters, while in reality they result from asymmetric 
mergers as well as from converging flows in large scale sheets and filaments. 
The more realistic emission from these complex geometries can be best modeled 
through detailed hydrodynamic simulations ]6l| , |6^ , ^6| . Mergers of comparable 
mass clusters, for example, would tend to produce only mild shocks due to the 
prior heating of the shocked gas, leading to reduced non-thermal emission due 
to the steep power-law slope of the accelerated electrons. 

In analogy with the collisionlcss shocks of supernova remnants, the inter- 
galactic shocks are also expected to accelerate a hadronic cosmic-ray component 
that would acquire a fraction (~ 10-50%) of the post-shock energy more sub- 
stantial than that carried by relativistic electrons (~ 1-10%, see JOI). Collisions 
of the hadronic cosmic-rays with IGM protons produce pions, w, which decay 
into 7-ray photons; however for typical parameters this radiation component 
is expected to be sub-dominant relative to the inverse-Compton production of 
7-rays by the relativistic electrons |^ , B(|p2) . The buoyancy of the relativistic 
cosmic-ray fluid may lead to instabilities that would eventually separate it from 
the IGM gas and allow it to expand into the large-scale voids that fill most of 
the volume in between the sheets and filaments of the IGM. For a relativistic 
fluid, the adiabatic decrease of the pressure with increasing volume (p oc V~ 4 / 3 ) 
is milder than for the IGM gas (p oc V~ 5 ^ 3 ). It is possible that the intergalactic 
voids are currently being dominated by cosmic-ray and magnetic field pressure 
which suppresses the formation of galaxies there. 
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